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ABSTRACT 

O 

We investigate energetic type Ic supernovae as production sites for 6 Li and 
Be in the early stages of the Milky Way Recent observations have revealed 
that some very metal-poor stars with [Fe/H]< —2.5 possess unexpectedly high 
abundances of 6 Li. Some also exbihit enhanced abundances of Be as well as N. 
From a theoretical point of view, recent studies of the evolution of metal-poor 
massive stars show that rotation-induced mixing can enrich the outer H and He 
layers with C, N, and O (CNO) elements, particularly N, and at the same time 
cause intense mass loss of these layers. Here we consider energetic supernova 
explosions occurring after the progeniter star has lost all but a small fraction of 
the He layer. The fastest portion of the supernova ejecta can interact directly 
with the circumstellar matter (CSM), both composed of He and CNO, and induce 
light element production through spallation and He-He fusion reactions. The 
CSM should be sufficiently thick to energetic particles so that the interactions 
terminate within its innermost regions. We calculate the resulting 6 Li/0 and 
9 Be/0 ratios in the ejecta+CSM material out of which the very metal-poor stars 
may form. We find that they are consistent with the observed values if the mass 
of the He layer remaining on the pre-explosion core is ~ 0.01 — 0.1 M , and 
the mass fraction of N mixed in the He layer is ~ 0.01. Further observations of 
6 Li, Be and N at low metallicity should provide critical tests of this production 
scenario. 
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Subject headings: nuclear reactions, nucleosynthesis, abundances — supernovae: 
general — stars: abundances 

1. Introduction 

Among the light elements Li, Be and B (LiBeB), 7 Li is thought to arise from a variety 
of processes, including big bang nucleosynthesis (Spite & Spite 1982), asymptotic giant 
branch stars, novae (D'Antona & Matteucci 1991), and the //-process in type II supernovae 
(Woosley et al. 1990); the latter may also contribute to n B. On the other hand, the main 
production channel for the rest, in particular for 6 Li and Be, is believed to be cosmic-ray 
induced nuclear reactions. The most widely discussed models of LiBeB production are based 
on cosmic rays accelerated in supernova shocks (Meneguzzi et al. 1971; Vangioni-Flam et 
al. 2000). Observations of metal-poor stars in the Galactic halo show a primary relation 
between [Fe/H] and [Be/H] or [B/H], which is consistent with spallation by cosmic rays 
enriched with C, N, and O (CNO) from fresh SN ejecta impinging on interstellar H or He 
(Duncan, Lambert, & Lemke 1992). 

An alternative possibility was proposed by Fields et al. (1996, 2002), who considered 
explosions of Type Ic supernovae (SNe Ic) as a site for primary LiBeB production. Since 
it is expected that a fraction of C and O in the surface layers of the ejecta are accelerated 
to energies above the threshold for spallation reactions when the shock passes through the 
stellar surface, LiBeB production can occur through the direct interaction of the ejecta with 
the ambient material, without the need for shock acceleration of cosmic rays. This was 
explored in greater depth by Nakamura & Shigeyama (2004), who used more realistic stellar 
models and equations of state, together with a 1-dimensional relativistic hydrodynamic code 
to investigate in detail how much of the ejecta mass acquires sufficient energies for the 
LiBeB production. All of these studies used stellar models which completely lost their H 
and He envelopes and assumed the target to be interstellar matter (ISM) consisting of 90% H 
and 10% He, ignoring any circumstellar matter (CSM). Therefore, the only reactions under 
consideration were C, O + H, He — > LiBeB. 

Recent observations by VLT/UVES (Asplund et al. 2005) and Subaru/HDS (Inoue 
et al. 2005, Aoki et al, in prep.) have revealed that some very metal-poor stars possess 
surprisingly high abundances of 6 Li, much higher than expected from standard supernova 
cosmic ray models (Ramaty et al. 2000; Suzuki & Yoshii 2001; Prantzos 2006). The measured 
values are also higher than can be accommodated in the SN Ic production scenario discussed 
above, even in the case of an energetic explosion similar to SN 1998bw, as can be seen from 
the results of Nakamura & Shigeyama (2004). However, besides spallation of C and O, the 
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fusion reaction He + He — > 6 Li may also be potentially important. This reaction will be 
significant in the conceivable case that a small fraction of He is still remaining in the surface 
layer of the core at the time of the explosion, while most of the He has been transported 
to the CSM through mass loss. Meynet & Maeder (2002) and Meynet et al. (2006) have 
recently calculated the evolution of metal-poor massive stars taking into account the effects 
of rotation-induced mixing, and shown that extensive mass loss occurs even in extremely 
metal-poor cases, in addition to significant enrichment of the N abundance near the surface. 
If this N is accelerated at the shock breakout of the SN explosion, a significant amount of Be 
can be produced through the reaction N + He — ► 9 Be, because of its low threshold E th (~ 6 
MeV/A) and high cross section at peak a p (~ 24 mb) compared with other 9 Be-producing 
reactions (eg. E th ~ 8 MeV/A and er p ~ 8 mb for O + He — > 9 Be, see Fig. 1). Indeed, 
recent observations indicate that the abundances of both Be and N may be enhanced in some 
metal-poor stars in which 6 Li is detected (Primas et al. 2000a,b; Israelian et al. 2004). In 
this Letter, we focus on the interactions He, CNO + He, N — > LiBeB occurring between the 
ejecta accelerated by the energetic SN explosion and the CSM, in order to account for the 
abundances in such very metal-poor stars. 

We summarize recent observations of the relevant elements in very metal-poor stars 
in §2. In §3 we discuss SN explosions of stars with stripped He layers embedded in the 
CSM. Our calculations and results are shown in §4 with an appropriate parameter range. 
Conclusions are presented in §5. 

2. Summary of Observations 

Observational determination of the abundances of 6 Li in metal-poor stars is extremely 
challenging due to its deficiency and the proximity of its absorption lines to the much stronger 
ones of the 7 Li isotope. Owing to the advent of large telescopes and the improved capabilities 
of observational instruments, reliable measurements have become possible only in the last 
several years (for earlier observations, see Smith et al. 1998; Hobbs et al. 1999; Cayrel et 
al. 1999). Asplund et al. (2005) observed 24 metal-poor halo stars with VLT/UVES and 
positively detected 6 Li in nine of them at the > 2a significance level. The subdwarf LP 
815-43 is the most metal-poor object ([Fe/H] = —2.74) with 6 Li in their sample. Inoue et al. 
(2005) reported a tentative detection with Subaru/HDS in the even more metal-poor star G 
64-12 ([Fe/H]= —3.17), although a detailed analysis is still ongoing (Aoki et al., in prep.). 
With VLT/UVES, Primas et al. (2000a,b) have measured the Be abundances of these two 
stars to be loge(Be) = -1.09 ± 0.20 for LP 815-43 and -1.10 ± 0.15 for G 64-12. The Be 
abundance in G 64-12 is considerably higher than that expected from previous measurements 
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of the Be-Fe relation in stars with similar metallicities, and this may be the case for LP 815- 
43 as well. Israelian et al. (2004) analyzed nitrogen abundances in 31 metal-poor stars and 
found that both LP 815-43 and G 64-12 are more N rich than average. The abundances of 
the relevant elements in these two stars are listed in Table 1. We adopt the values obtained 
by 1-D LTE analyses to ensure consistency. 

3. Metal-Poor SN Ic Embedded in CSM 

Meynet et al. (2006) simulated the evolution of metal-poor ([Fe/H] =— 6.6, —3.6) massive 
(60 M Q ) stars with rapid rotation (800 km/s) and found that C is enhanced in the outer 
layers by rotation-induced mixing, promoting intense stellar winds and significant loss of 
their envelopes in spite of their initial low metallicities. At the same time, the mixing results 
in enrichment of N in the He layers. If a fraction of the He layer is still remaining when 
the SN explodes, it is expected that the N within is accelerated at the shock breakout, and 
a significant amount of Be will be produced through the reaction N + He — ► 9 Be, thanks 
to its low threshold and high cross section at peak. Thus, we consider stars that have 
lost all of their H layer and most of their He/N layer before the explosion. Here we use 
an explosion model of a ~ 15 M Q core, originating from a main-sequence star with mass 
Mnjs ~ 40 M & (Nakamura et al. 2001). The explosion energy is assumed to be 3 x 10 52 
ergs, corresponding to an energetic explosion similar to SN 1998bw. The mass of the ejecta 
becomes 13 M , containing 10 M of oxygen. The accelerated ejecta consisting of He and 
CNO will collide with the circumstellar He and N stripped from the progenitor star and 
produce LiBeB through the reactions He, CNO + He, N — > LiBeB. The energy distribution 
of the C/O ejecta before interaction is that calculated in Nakamura & Shigeyama (2004), 
shown in Figure 1 together with the cross sections for selected reactions (Read & Viola 1984; 
Mercer et al. 2001). Rather than recalculating the explosion hydrodynamics with the added 
He/N layer, we approximate by changing the composition of the accelerated outermost ejecta 
from C/O to He/N. This will not lead to considerable errors as long as the replaced mass is 
small. 

The "thick target" approximation is used, that is, the circumstellar He is so thick that 
light element production occurs entirely within the CSM while the ejecta particles lose energy 
mainly by Coulomb collisions with free electrons. This assumption is valid when the mass 
loss rate M is greater than 10~ 6 M yr _1 for the typical wind velocity of v w ~1,000 km 
s -1 , as can be seen by comparing the windblown material's mass column density a with the 
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stopping range R. We obtain for an a-particle with initial energy e 
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^ ~ 0.096 ( vlQMeV / A J ^i O -6M /yrJ \1, OOOkms" 1 ) \R & ) ' ^ 

which is well below unity. Here R is defined as the column depth through which a particle 
loses all of its energy. Even if the matter in the wind is not ionized, the energy loss rates 
due to ionization will be similar. Other processes such as escape from the system can be 
ignored. The LiBeB yields are calculated using the cross sections given by Read & Viola 
(1984) and Mercer et al. (2001). We also assume that the composition of the innermost 
CSM (i.e. the last wind material) concerned with the reactions is the same as that of the 
progenitor's outermost layers. 

As a result of the localized production of the light elements, their abundance ratios with 
respect to heavy elements averaged over the CSM and SN ejecta are likely to be inherited 
by stars of the next generation, as pointed out by Shigeyama & Tsujimoto (1998). Thus in 
the next section we compare the abundance ratios calculated from the above model with the 
metal-poor stars discussed in the preceding section. 



4. Origin of Li and Be in Metal-poor Stars 

We focus on the abundances of the very metal-poor star LP 815-43 from which 6 Li, Be 
and N have been detected, all at apparently enhanced levels. As seen in Table 1, Xe^/Xo ~ 
6.88l|^ x 10~ 7 and Xg Bc /Xo ~ 1-32+°;^ x 10~ 8 for this star. The mass M He , N of the He/N 
layer on the pre-explosion core and its mass fraction X N of N are the main parameters. Figure 
2 shows our results. The yields of 6 Li increase until M He ,N ~ 0.01 M , which corresponds 
to the threshold energy of the He + He — > 6 Li reaction (~ 11 MeV/A, see Fig. 1), and 
saturates for larger Mn e , n- They depend only slightly (~ 0.4%) on Xn (as long as Xn is 
not very large), because most of the 6 Li is produced through the He + He reaction. Only 
the line corresponding to the case of X^ = 0.005 is shown in the top panel of Figure 2. On 
the other hand, the yields of 9 Be, which is mainly produced through the C, O + He — > 9 Be 
reaction when the He/N layer is deficient, depend strongly on X^. Without N (Xn = 0), 
the Be yield decreases rapidly with M He ,N, since the He + He reaction does not produce Be. 
For small A N , Be yields first decrease with M He) N for the above reason, and then turn to 
increase since the cross section of the N + He — > 9 Be reaction peaks around ~ 13 MeV/A, 
corresponding to ~ 0.013 M integrated from outside for SN 1998bw model (see Fig. 1). 
It should be noted that if the He layer is sufficiently large, for instance Mh c > 1 M , our 
results here tend to overestimate the light element yields because a core with such a large 
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He layer might not be so compact in reality and effective acceleration at the shock breakout 
is not expected. To specify how much He can be accelerated above the production threshold 
in such cases, detailed calculations for the structures of mass-losing stars are required. 

Both Xe Li /Xo and X9 Bc /X show good agreement with the observational data of LP 
815-43 when M He ,N ~ 0.01 - 0.1 M and X N ~ 0.005 - 0.01. These are consistent with 
recent simulations of the evolution of metal-poor massive stars with rotation: Xn ~ 0.008 
for M ms = 40 M Q (Hirschi et al., in preparation) and ~ 0.01 for M ms = 60 M & (Meynet 
et al. 2006). The observed X.9 Bc /Xq for G64-12 can also be reproduced with parameters 
analogous to those for LP 815-43. The isotopic ratio 6 Li/ 7 Li in our models lies between 0.27 
and 1.4, which is much larger than the observed value 0.046 ± 0.022 for LP 815-43 (Asplund 
et al. 2005) so that the scenario discussed here is not expected to contribute significantly to 
7 Li. 



5. Conclusions and Discussion 

In this Letter we have proposed a new mechanism to produce the light elements, espe- 
cially 6 Li and Be recently observed in metal-poor stars. This is based on recent theoretical 
findings that rotating, metal-poor massive stars can lose substantial amounts of their en- 
velopes, and end up with SNe Ic with a small amount of He and N in the outermost layers. 
The outer layers composed of He and CNO are accelerated at the shock breakout and un- 
dergo spallation and He- He fusion reactions to produce mainly 6 Li and 9 Be. Our calculations 
show that if ~ 0.01 — 0.1 M of the He layer is remaining on the core in a SN explosion of 
a massive star with M ms > 30 M , then the observed abundance ratios 6 Li/0 can be repro- 
duced. If the He layer contains a small amount (~ 0.5 — 1%) of N, we can also reproduce 
9 Be/0. 

Compared with the standard picture of cosmic ray shock acceleration by normal super- 
novae, our high yield of 6 Li results from three crucial differences. First, we consider a large 
explosion energy, of order 10 52 erg, as observed in some SNe such as SN 1998bw. Second, 
our mechanism considers the direct interaction of fast ejecta with the CSM, and does not 
involve an efficiency factor for cosmic ray acceleration and is not affected by losses or escape 
during ISM propagation. Third, the energy distribution of the accelerated particles is a very 
steeply falling power-law, with spectral index ~ 4.6 as opposed to 2 for shock accelerated 
particles. This means that most of the energy is contained in the lowest energy portions, 
where the cross sections for both He + He — > 6 Li and N + He — > 9 Be peak (see Fig. 1). In 
fact, the total energy of particles above 10 MeV/nucleon in our fiducial model is 8 x 10 50 
erg, much higher than is achievable in the standard picture, 10-30% of 10 51 erg. 
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One may claim that 6 Li, and probably also 9 Be (Garcia-Perez & Primas 2005), in metal- 
poor stars may have been depleted from their initial values, since the observed abundances 
of 7 Li/H are a few times smaller than that predicted by standard big-bang nucleosynthesis, 
and 6 Li is more fragile than 7 Li (Asplund et al. 2005). The actual survival fraction of 6 Li 
is difficult to evaluate reliably because the pre-main-sequence Li destruction is sensitive to 
convection, which cannot be modeled without free parameters. However, depletion factors as 
large as 10 may still be compatible with our picture. The estimates in the preceding section 
is based on an explosion energy of the order of 10 52 ergs. The mass of ejecta with energy 
per nucleon above a certain value scales very strongly with the explosion energy, E cx as 
(Nakamura & Shigeyama 2004). Thus a SN with the explosion energy 2 times higher than 
that considered in the preceding section will be able to produce ~ 10 times larger amounts of 
light elements depending on the distribution of He and N. Nevertheless, if such very energetic 
explosions turn out to be rare, the mechanism proposed here may play a rather limited role 
compared to other potential LiBeB production processes. 

Both the degree of mass loss and the amount of N enrichment in the progenitor star 
are expected to be sensitive to its initial rotation speed (Meynet & Maeder 2002; Meynet 
et al. 2006). The actual rotation speed is presumably distributed over a wide range, as is 
the explosion energy and the extent of mass loss at the time of the explosion. Therefore, 
dispersions in the 6 Li and Be abundances are expected, and current observations suggest that 
this may indeed be the case for 6 Li (Aoki et al. 2004; Asplund et al. 2005; Inoue et al. 2005) 
as well as Be (Boesgaard & Novicki 2005). Our scenario predicts a close relation between Be 
and N, since Be arises directly as a consequence of N spallation. A looser correlation between 
6 Li and N is also expected, as effective 6 Li production requires sufficient mass loss, which 
in turn implies significant N enrichment. Further observations of 6 Li, Be and N for a large 
sample of metal-poor stars should provide definitive tests. Note that such correlations are 
not expected for other scenarios which involve mainly 6 Li production only, such as structure 
formation cosmic rays (Suzuki & Inoue 2002; Rollinde et al. 2005), active galactic nuclei 
outflows (Prantzos 2006; Nath et al. 2006), and production processes in the early universe 
(Jedamzik et al. 2005). It is also mentioned that mass loss onto companion stars in binary 
systems may represent an additional pathway to our scenario, provided that sufficiently thick 
CSM is remaining at the time of the explosion. We reiterate that more observational data is 
necessary to elucidate what fraction of the LiBeB abundances seen in halo stars of different 
metallicity can be explained by the mechanism proposed here. 

We are grateful to Georges Meynet, Raphael Hirschi, Takeru K. Suzuki and Sean Ryan 
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which led to clarification of a number of points in our original manuscript. This work has been 
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Table 1. Abundances from the literature. 



Star 


[Fe/H] 


loge( 7 Li) 


6 Li/ 7 Li loge(Bc) 


loge(N) 


loge(O) 


LP 815-43 


-2.74 a 


2.16 b 


0.046±0.022 b -1.09±0.20 c 


5.61 a 


6.54 c 


G 64-12 


-3.17 d 


2.30° 


-1.10±0.15 c 


6.15 a 


6.34 d 



a Israelian et al. (2004) 
b Asplund et al. (2005) 
c Primas ct al. (2000a) 
d Akerman ct al. (2004) 
e Primas ct al. (2000b) 
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Fig. 1. — The energy distribution of ejecta and cross sections of reactions as functions 
of energy per nucleon e. The solid line represents the result of a numerical calculation 
(Nakamura & Shigeyama 2004) where M(> e) denotes the mass of ejecta that have particle 
energy per nucleon greater than e. The dotted, dashed, and dash-dotted lines show the cross 
sections a of reactions He + He — > 6 Li, N + He — > 9 Be, and O + He — > 9 Be, respectively. 
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Fig. 2. — Mass ratios of 6 Li (top panel) and 9 Be (bottom) to O as functions of M HCj n- 
Shown are the cases with mass fractions of N X N = 0.005 (solid line), X N = 0.01 (dotted), 
X N = 0.02 (dash-dotted), and without N (X N = 0; dashed). Shaded regions represent the 
observed ratios for the very metal-poor star LP 815-43 (Primas et al. 2000a; Asplund et al. 
2005) including the error. 



